The Nature of the X-ray Source in SNR RCW 103 
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ABSTRACT 

The discovery of the 6.67 hr periodicity in the X-ray source IE 161348—5055 
associated with the supernova remnant RCW 103 has raised interesting sugges- 
tions about the nature of the X-ray source. Here we argue that in either accreting 
neutron star or magnetar model, a supernova fallback disk may be a critical in- 
gredient in theoretical interpretations of IE 161348—5055. We further emphasize 
the effect of fallback disks on the evolution of young compact objects in vari- 
ous ways, and suggest that even SS 433 could also be powered by fallback disk 
accretion process. 
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Subject headings: stars: individual (IE 161348—5055, SS 433) - stars: neutron - 
supernova remnants - supernovae: individual (RCW 103)- X-rays: stars 



1. Introduction 

The enigmatic radio-quiet central compact objects discovered in supernova remnants 

(SNRs) have challenged conventional tho ughts that most young neutron stars (NSs) evolve in 

a ma nner similar to Crab-like pulsars (see IPavlov. Sanwal. fc Teterkooi iKaspi. Roberts. Sz Harding! 
20061 . for recent reviews). The X-ray source IE 161348—5055 (hereinafter 1E1613), the pro- 
totype of the growing radio-quiet objects in SNRs, was first detected with Einstein as a 



faint , unresolved source located n ear the center of the SNR RCW 103 (ITuohv fc Garmire 
1980 ) with an age of ~ 2000 yr (INugentl Il984l) and distance of ~ 3.3 kpc ( jCaswell et al. 



19751 ). From X-ray observations with ASCA, iGotthelf. Petre. fc Hwang! (119971 ) found a com 



pact X-ray source inside RCW 103 with the X-ray luminosity Lx ~ 10 34 ergs -1 and the 
black-body temperature about 0.6 keV. Further ASCA observation of 1E1613, together with 
archived Einstein and ROSAT data showed t hat this source manifested an o rder-of mag- 
nitude decrease in luminosity over four years (IGotthelf. Petre. fc Vasishtl Il999l ) . suggesting 
that this object may be an accreting source rather a cooling NS. But no X-ray pulsations, 
radio or optical counter part were detected. From Chandra observations and archival ASCA 
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data, I Gar mire et al.l (120001 ) reported the ~ 6 hr periodicity of its flux. The lightcurve taken 
by XMM-Newton did not sh ow the ~ 6 hr periodicity but discovered an eclipse about 3 hr 
( IBecker fc Aschenbachll2002l ). More recent observatio ns with XMM-Newto n showed a strong, 
unambiguous periodic modulation at 6.67 ±0.03 hr (IDe Luca et al.ll2006l ). 



The origin of this period could be either an orbital or a spin period. In the former 
case 1E1613 would be a low-mass X-ray binary (LMXB) which survived the SN event; 
i n the latter, 1E161 3 could be an isolated magnetar with magnetic field of order 10 15 G 
( ]De Luca et al.ll2006l ). In each possibility the source remains to be peculiar and puzzling. In 
this Letter, we discuss the nature of the X-ray source 1E1613 and its implications for young 
compact objects, arguing that in both scenarios a fallback disk may be required to account 
for the observational characteristics of 1E1613. We further propose that fallback disks may 
play a critical role in the evolution of several kinds of young compact stars including the 
famous X-ray source SS 433. Throughout this study we assume that both 1E1613 and RCW 
103 originated from the same SN event. 



2. Accreting NS interpretation 



The long term flux change suggests that 1E1613 may be an accreting source (IGotthelf. Petre. &: Vasisht 



19991 ). The lightcurve observed by XMM-Newt on shows for the evidence o f an eclipse, in- 
dicating IE 1613 as a NS in a binary system (IBecker fc Aschenbachl |2002| ). The possible 
companion should be less massive than 0.4 M^, if it is a normal star, as indicated by op - 
tical/IR observations (jPavlov. Sanwal. fc Teterl 12004 ; IvVang. Kaplan, fc Chakrabartyl 120071 ). 
Thus 1E1613 could be a young NS accreting from a low-mass companion star. 

The LMXB interpretation meets many difficulties. One of them is the very low birthrate 
of LMXBs. It is well known that LMXBs are extremely rare objects, with a birthrate of 
< 10~ 6 — 10~ 7 yr _1 . The birthrate of 1E1613 as an LMXB would be even lower due to the very 
low initial mass of its companion - the proge nitor binary is very likely to either merge dur - 
ing the previous common envelope evoluti on (IJustham. Rappaport. fc Podsiadlowskill2006l ). 
or be disrupted during the SN explosion ( Brandt fc Podsiadlowski 19951 ). This makes the 
probability to find an LMXB like 1E1613 inside a SNR extremely low. Moreover, if 1E1613 
is indeed an LMXB, it is difficult to explain why there is few association between SNRs and 
intermediate- or high-mass X-ray binaries (HMXBs) in the Galaxy, which are more likely to 
form and to be observed. Only a single SNR-HMXB association (SS433 and the SNR W50, 
Geldzahler. Pauls, fc Salter! Il980l ) is known in the Galaxy (see discussion below). 



Even if the LMXB can survive the common envelope evolution and the SN explosion, 
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it will take at least ~ 10 8 — 10 9 yr for orbital decay so that Roche-lobe overflow from 
the companion starts (the companion star of mass ~ 0.4 M will take a time longer than 
the age of the Galaxy by nuclear evolution to fill the Roche-lobe). This timescale is in 

< 2000 yr) of the SNR RCW 103 



sharp contrast with the young age (~ 2000 yr) of the SNR RCW 103. iDe Luca et al. 
( 120061 ) suggest that the LMXB may have a significant orbital eccentricity (due to the SN 
kick) and hence orbital modulation in the captured mass by the NS from the companion's 
wind. We note that the stellar wind mass loss rate for a < 0.4 M Q dwarf is as low as 
~ 10~ 15 M yr -1 , and obviously cannot power the X-ray luminosity (~ 10 33 — 10 35 ergs" 1 ) of 
1E1613. Thus the wind from the donor star, if really relevant to the flux modulation, should 
be induced by X-ray irradiation from the NS. X-ray irradiation can affect stars of mass 
< 1.5 M Q under certain circumstances by ionizing the hydr ogen at the base of t he irradiated 
su rface layer and disrupting the surfa c e con vection zone (IPodsiadlowskil Il99ll ). According 
to lPfahl. Rappaport. fc Podsiadlowskil (120031 ). the maximum orbital period for which X-ray 



irradiation is important is 



orb, max 



70days(^W 



S c ,n 

is the accretion rate onto the NS in units of 10~ 8 1 1 1, -1 



Mq yr i , S C) n is the critical X- 
2 , and e < 1 is a factor that takes 



where M. 

ray flux for hydrogen ionization in units of 10 11 ergs _1 cm 
into account the geometry of the accretion disk and the star, th e albedo of the star, a nd the 
fraction of X-rays that penetrate below the stellar photosphere ( iHameury et al.lll993l ). Take 
L x ~ 10 34 ergs -1 for 1E1613, Eq. (1) yields P rb,max < 1-7 hr, still smaller than the 6.67 hr 
period. However, if the orbit is highly eccentric, intensive winds might be excited by X-ray 
irradiation at periastron. Nevertheless, if future observations confirm the existence of the 
binary companion, it will have important impact on the current theories of NS formation. 

In our opinion, perhaps it is more natural to assume that 1E1613 is a NS accreting from 
the fallback disk rather a low-mass companion. Current stellar evolution models predict that 
during the core collapse of massive stars, a considerable amount of the stellar material will fall 
back ont o the compact, collapsed remnants (NSs or BHs), usually in the form of an accretion 
disk (e.g. IWoosleylll993l ). This point of view is sup ported by the discovery of a remnant disk 
around the anomalous X-ray pulsar 4U 0142+61 (jWang. Chakrabarty. fc Kaplan! 120061 ). In 
this case one does not need to care about the age contrast in the low-mass donor and 
young SNR in the LMXB scenario. The lifetime of the fallback disks is around 10 4 — 
10 5 yr. Du ring this time the mass supply may be s ufficient to power the observed X-ray 



luminosity (IChattterjee. Hernquist. &: Narayanl 120001 ). The flux outbursts and dips in the 
l ight curve could be d ue to disk instabilities and occultations by disc structures, respectively 
(IDe Luca et al.ll2006l). The 6.67 hr perio dicity might be rel ated to the precession of the NS 
with a hot spot (IHeyl fc Hernquistl 120021 ) or the disk itself (jKatzlll973l ). 
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Popovl (120071 ) suggested an interesting idea that 1E1613 could be in a double NS binary 
system, where the newborn NS which produced the SNR has a remnant disk around it, and 
the older NS, identified as 1E1613, accretes periodically from that disk when passing close 
to the companion. The situation here is similar to what happens in Be/X-ray binaries, in 
which a NS accretes from the disk-li ke winds from its Be companion star in an ecce ntric 
orbit. In the disk truncation model (jOkazaki fc Negueruelall200ll ; lOkazaki et al.l 120021 ) the 
X-ray outbursts in Be/X-ray binaries are explained as follows. The NS exerts a negative tidal 
torque on the viscous decretion disk of the Be star, resulting in the truncation of the disk. 
The disk matter would then accumulate in the outer rings of the disk until the truncation was 
overcome by the effects of global one-armed oscillations, disk warping, etc. The subsequent 
sudden infall of high-density disk matter onto the NS causes type II X-ray outbursts. If the 
tidal truncation is not very efficient and the disk extends beyond the Roche-lobe of the Be 
star at periastron, the matter could be accreted onto the N S during the periast r on pa ssage, 
resulting in (quasi-)periodic type I bursts. Calculations by IZhang. Li. fc Wang! (120041 ) have 
shown that the narrower the binary system, the more efficient the truncation, since the 
truncation efficiency oc P^b ■ The orbital period (6.67 hr) of 1E1613, which is much smaller 
than those of Be/X-ray binaries, suggests that the flux variations in 1E1613 may be erratic 
rather periodic, if caused by the disk truncation effect. 

One of the problems for all the accreting NS scenarios is that, to produce an accretion 
luminosity in the range observed, the NS would have a low magnetic field and/or a slow 
rotation period, i.e., P ~ (0.35 — 2.5)£>^q 7 s, where B 10 is the NS magnetic field in units of 
10 10 G, so that the accreting material can overcome the centrifugal barrier (IDe Luca et al. 



20061 ). Although these values seem to be peculiar if compared with the canonical values 
(B ~ 10 12 G and P ~ 10 — 100 ms) of young NSs, there is accumulati ng; evidence for some 
NSs born spinning slowly and with a relatively weak magnetic field (IHalpern et al.l 120071 : 
Gotthelf fc Halpernll2007l l 



3. Magnetar interpretation 



De Luca et al.l (120061 ) favor the idea that 1E1613 may be a magnetar rotating at 6.67 



hour with B > 10 15 G, as the X-ray variabilities, luminosity and spectral shape would be 
naturally explained in the magnetar frame. Even with such field a NS cannot spin down to 
6.7 hours via magneto-dipole radiation (MDR) during the lifetime of the SNR. Therefore, a 
SN fallback disk is requi red to interact with th e NS, providing additional propeller spin-down 
torque. Calculations by IDe Luca et al.l (120061 ) show that a disk of 3 x 10~ 5 M would have 
been enough to slow down, over 2000 years, a5 = 5x 10 15 G magnetar to current period, 
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provided that the initial spin period should be Pi > 300 ms. Magnetic field in magnetars 
is generally thought to be generated by turbule nt dynamo, whose strength depends on the 
star's rotation rate (IDuncan fc Thompson! Il992l . see however, Vink & Kuiper 2006). Such a 
period seems to be too long for magnetars, although its possibility cannot be ruled out. 

Additionally, the known magnetars include the anomalous X-ray pulsars (AXPs) and 
s oft g amma-ray repeaters (SGRs), which are rotating at ~ 5 — 10 s (IWoods fc Thompson 



20061 ). If 1E1613 is a magnetar, one needs to explain why its period is much longer than 
those of AXPs and SGRs. 

The mode of interaction between a NS and a surrounding disk is determined by the 
location of the inner radius R in of the disk with respect to the characteristic radii, the 
corotation radius R c = (G MP 2 / 'An 2 ) 1 / 3 , and the light c ylinder radius = cP/2tt, where 



M is the mass of the NS (e.g. Illlarionov fc Sunyaevill975l ; lLipunovlll992l ). The position of the 



inner radius of the disk can be estimated by comparing the electromagnetic energy density 
generated by the NS with the kinetic energy density of the disk. The NS is expected to be 
in the propeller and ejector (radio pulsar) stage if the inner radius of the disk is beyond the 
corotation and light cylinder radius, respectively. Since the kinetic energy density in the disk 
has the dependence oc r -5 / 2 on the radial distance r from the center of the NS, steeper than 
the electromagnetic energy density (or radiation pressure) outside the light cylinder (oc r -2 ), 
stable equilibrium of the disk ou tside the light cylinder is not allowed, unless it is beyond 
the gravitational capture radius (lLipunovlll9921 ). 



In their calculations iDe Luca et al.l (120061 ) have adopted the traditional estimates of the 
radiation pressure from a NS by em ploying a rotat i ng ma gnetic dipole in vacuum to generate 
the electromagnetic fields. Recently lEksi fc Alparl (120051 ) derived t he ele ctromagnetic energy 
density from the global electromagnetic field solution of iDeutschl (119551 ) for obliquely rotat- 
ing magnetic dipoles. They showed that the electromagnetic energy density of a rotating 
dipole makes a rather broad transition for disk existence across the light cylinder for small 
inclination angles. Within this model the disk could survive beyond the light cylinder even if 
the radio pulsar activity turns on. 

To examine to what extent the fallback disks affect the spin evolution of magnetars, 
we carried out Mon te Carlo s imula tions of the evolution of 10 6 NSs based on the spin-down 
model presented in iLipunovl (119921 . see also Li 2002). As the mass of a fallback disk is not 
replenished, mass flow rate in the disk declines and the inner radius of the disk moves out. 
The NS in this case generally passes three evolutionary stages. (1) First is the "ejector" 
phase, in which the radiative pressure from the NS is sufficient to keep the surrounding 
plasma away from the light cylinder; the NS evolves as a radio pulsar. Here we assume that 
once i?i n /i?L > Xct, the gaseous disk is disrupted by the radiation pressure, and the NS is 
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spun down only by MDR (her e x„ > 1 is a param eter depending on the magnetic inclination 



and is calculated according to lEksi &: Alparll2005l ); if 1 < R^/R^ < x cr , the disk is still stable 



beyond R^, and MDR also works. (2) If R c < R in < R L , the "propeller" phase follows, in 
which the plasma interacts with the neutron star magnetosphere but further accretion is 
inhibited by the centrifugal barrier, and the disk exerts a propeller spin down torque on the 
NS. (3) For sufficiently long time of evolution, the NS will enter the "accretor" phase, in 
which the inner disk radius becomes smaller than the corotation radius and steady accretion 
of the plasma is allowed. 

In our calculations the initial NS magnetic field B is chosen from a log normal distri- 
bution of mean 15 and standard deviation 0.4. We assume that all NSs were born with a 
surrounding supernova fallback disk, the initial masses of which log(AM/M ) are distributed 

uniformly between —6 and —2. Mass flow rate through the disk is assumed to decline in a 

powe r law with time, M oc r L25 JCannizzo. Lee. & Goodm^r]|l99ol : iFrancischelli. Wijers. fc Brown! 



20021 ). We set the initial spin periods Pi and the inclination angle to be distributed uniformly 
between 2 ms and 50 ms, and between 0° and 90°, respectively. For the propeller torque, we 
adopt the expression T prop = MRf n [Qk (-Rin) — &s), where Ok (- Rin) is the Kepleria n angular 



J. J. I"- Ill L \ / " J ' \ ■ I I I |i . . I . «_< 

velocity at -R in and Vt s the angular velocity of the NS (see also iMenou et al.lll999l ). We stop 



the calculations at a fiducial time of 2500 yr, to be compatible with the age of RCW 103. 

In Fig. 1 we plot the histogram of the spin periods for the magnetars. The hatched 
regions indicate the propeller/accretor (+45°) and ejector (—45°) systems, respectively. It is 
noted that most (~ 99%)Q of the magnetars are in the ejector state, having spin periods of a 
few seconds. These features are consistent with most of AXPs and SGRs, which occupy the 
majority of the magnetar family. The spin periods of the propeller/accretors, which consist 
of ~ 0.6% of the whole population, have a much broader distribution, peaking around a 
few of 10 3 s. 1E1613 may belong to this latter group. In this case, the puzzling initial spin 
period P > 300 ms does not appear. 



4. Discussion 

The above arguments suggest that in either accreting NS or magnetar interpretation, a 
SN fallback disk seems to be favored, to explain the observational characteristics of 1E1613. 
Especially we show that, in the more reasonable magnetar model, it is possible to model 
the spin period evolution in 1E1613 with typical values of the input parameters for a mag- 



lr The percentages presented here are subject to the uncertainties in the initial parameters and the spin- 
down mechanisms adopted, and should not be taken very seriously. 
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netar. This work emphasizes that fallback disks may considerably influence the formation 
and evolution of young compact objects evolved from core collapse of massive stars. The 
existence of fallb ack disks can be tested by searching for optical/IR emission from the coole r 
parts of t he disk (IPerna. Hernquist. &: Naravanll2000l ; IWang. Chakrabarty. &: Kaplan! 120061 ). 
Recently IWang. Kaplan, fc Chakrabarty! (120071 ) reported on search for the optical/infrared 
counterparts to the central compact objects in four young SNRs including RCW 103, but 
found that there is confusion with several faint stars at the position of 1E1613. 

The fallback disk may manifest itself in various ways. Recent X-ray observations show 
that some young pulsars, such as the Crab and Vela pulsars, may have the jet configuration , 



which suggests the existence of a disk surrounding the neutron star (IBlackman fc Pernall2004l . 
and references therei n) . Such disks can influence the braking indices and timing ages of some 
young radio pulsars (IMenou. Perna. fc Hernquistll200ll ; iMarsden. Lingenfelter. fc Rothschild 
2001 ; Jiang fc Li 2005h . A sm a ll gro up of rotating radio transients (RRATs) were recently 
reported by iMcLaughlin et al.l (120061 ). These objects are charac terized by single, dispersed 
bursts of radio emission with durations between 2 and 30 ms. |LJ (120061 ) suggested that these 
phenomena could be due to the interaction between the NS magnetosphere and the sur- 
rounding SN debris disk. The pattern of radio pulsar emission depends on whether the disk 
can penetrate the light cylinder and efficiently quench the processes of particle production 
and acceleration inside the magnetospheric gap; a precessing disk m a y naturally account 
for th e switch-on/off behavior in PSR B1931+24 (IKramer et al.ll2006l ). ICordes fc Shannon 
( 120061 ) also proposed that the bursty emission by RRATs is due to a disk of circumpulsar 
asteroids randomly straying into the magnetosphere. 

Fallback disk accretion may also play an important role in producing the X-ray emis- 
sion of some ultral uminous X-ray sourc es (ULXs) that are possibly associated with SNRs, 
e.g. IC 342 X-l ([Roberts et all l2003h and the ULX in the SNR MF 16 in NGC 6946 
([Roberts fc Colbertl 120031 ). This association seems to rule out mass transfer in binaries 
as the main energy source, since there is little time for the donor stars (if they exist) to 
evolve and transfer mass rapidly to the black holes. A much more plausible explanation is 
tha t the black hole is accreting from the disk originating from the fallback supernova debris 
(H l2003h . 



In our Galaxy the X-ray source SS 433 is usually thought to appear as an ULX if 



observec 


face- 


Mareon 


1984; 



the presence of its relativistic jets precessing with a periodicity of about 162 d. There is no 
consensus about the binary component masses and the origin of W5 0. It might be a stellar 



windblown bubble produced by the SS 433 jets, as was proposed by iBegelman et al.l (11980 
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see however, Zealy et al. 1980). The companion sta r has been suggested to be a supergiant 
which is filling its Roche lobe (jHillwig et al.l 12004 . but see Barnes et al. 2006), thereby 
producing an extre mely high-mass transfer that occu rs on thermal time scale with the rate of 
M ~ 10~ 4 Mq yr _1 (IKing. Taam. &: Begelmanll2000l ). This gives rise to a total luminosity of 

10 35 -10 36 



10 40 ergs 1 , although the observed unabsorbed 2 — 10 keV X-ray lum inosity is 



ergs l . However, recent Chandra observations by lLopez et al.l (120061 ) suggest the radius of 



the compa nion to be about one third of the Roche lobe radius. A similar conclusion was also 
reached by iBrinkmann. Kawai. fc Matsuokal (119891 ) with Ginga observations. If it is correct, 
this indicates that the compact object (possibly a black hole) is likely to be accreting through 
a SN fallback disk rather Roche lobe overflow. Stellar wind accretion is also possible. But it 
is difficult to account for the high mass accretion rate and the existence of a stable accretion 
disk. Additional evidence for fallback disk accr e tion i n SS 433 my lie in recent XMM-Newton 
observations by IBrinkmann. Kotani. fc Kawail (120051 ). These authors find over- abundances 
in different elements in the outflowing gas: Si and S by factors of ~ 2, Ni by ~ 8. This 
overabundance of heavy elements might be a tracer for the past explosion of the massive 
progenitor of the compact star, since elemental synthesis of for example Ni in the inner parts 
of the accretion disk appears unlikely due to the extremely high temperatures required. 
The fallback disk accretion is the most natural way to convert the SN ejecta containing the 
synthesized heavy elements into the outgoing jets. 

In summary, SN fallback disks may account for the observed features of young compact 
stars through various ways of interaction between the disk and the star. Our results motivate 
further efforts to detect disks around NSs and BHs in SNRs, and more detailed models of 
the disk-star interaction. 



This work was supported by the Natural Science Foundation of China under grant 
numbers 10573010 and 10221001. 
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Fig. 1. — Distribution of the spin periods of magnetars at the age of 2500 yr. The hatched 
regions indicate propeller/accretor (+45°) and ejector (—45°) systems, respectively. 



